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ABSTRACT 

The study of dynamic stall of rapidly pitching airfoils is being 
conducted at NASA Ames Research Center. Understanding this 
physical phenomenon will aid in improving die maneuverability 
of fighter aircraft as well as civilian aircraft A wind tunnel 
device which can linearly pitch and control an airfoil with rapid 
dynamic response is needed for such tests. 

To develop a mechanism capable of high accelerations, an accu- 
rate model and control sy stem is created. The model contains 
mathematical representations of the mechanical system, includ- 
ing mass, spring, and damping characteristics for each structural 
element, as well as coulomb friction and servovalve saturation. 
Electrical components, both digital and analog, linear and 
nonlinear, are simulated. The implementation of such a high- 
performance system requires detailed control design as well as 
state-of-the-art components. 

This paper describes the system model, states the system 
requirements, and presents results of its theoretical performance 
which maximizes the structural and hydraulic aspects of this 
system. 

INTRODUCTION 

Assorted studies have investigated dynamic stall and its relation 
to rapidly pitching airfoils. Results indicate that pitch angles sig- 
nificantly beyond the static stall angle can be reached with large 
excursions in lift and pitching moment. The large excursions in 
lift and pitching moment are, in fact, the primary limiting factors 
in the performance of many aircraft (1). 

Previous airfoil pitching mechanisms have achieved pitch rates 
of 800 deg/sec at an average power of 0.6 hp. The required 
dynamics of the present experiment are 3600 deg/sec at an aver- 
age power of 10 hp with peak power exceeding 40 hp. To per- 
form the required tests with an airfoil, the pitching mechanism 
must meet the following maximum dynamic parameters: 

pitch angle, 0-60° 
velocity, 3600°/sec 
acceleration, 600,000°/sec 2 
acceleration time, 0.006 sec 
acceleration angle, 10° 

The large acceleration required at Mach 0.5 airflow necessitates a 
powerful prime mover, or energy source. Detailed discussion of 
the prime mover selection will be omitted in this paper, however, 
an agile, high-performance hydraulic system provided the best 
power and response for its physical size. The emphasis of this 


discussion will be on control-system modeling and design. This 
discussion will include control theory and mathematical device 
modeling, followed by appropriate Bode and root locus plots to 
demonstrate system performance. The design was accomplished 
using the control system analysis program MATRIXx/System 
Build (Integrated Systems Incorporated). Each of the system 
components will be represented by a mathematical model, or 
transfer function, expressed in the complex frequency domain. 
Time-dependent solutions will be presented in the form of 
response plots. 

CONSTANT PITCH RATE ACTUATOR 

The actuator is shown in Figure 1 installed on the wind tunnel 
test section. The airfoil is mounted between two glass disks 
which rotate about the quarter cord of the airfoil. Glass disks 
provide optical access for measuring the vortices developing 
across the airfoil's surface. 

The actuator is composed of a hydraulic servovalve, cylinder, 
and truss structure. A hydraulic servovalve operates a cylinder 
which pulls on a truss structure. The truss reacts on a bellcrank 
which pulls connecting links vertically, thereby pitching the 
glass disks and the airfoil. The torsion tube in the bellcrank is 
needed to minimize the torsional moment across the small air- 
foil. The NACA 0012 profile airfoil to be used is small 
(9.8 x 3.0 in.); therefore, the slightest twisting moment 
across the airfoil will disrupt the developing vortices. 

SYSTEM MODELING 

To study the feasibility and eventually maximize the pitching 
mechanism, a model representing the physical components is 
created. This model includes electrical components such as the 
transducers, motion controllers (digital and analog), and the 
software routines contained within the digital controller. The 
mechanical structure is modeled to include mass, spring, damp- 
ing characteristics, and coulomb friction. Finally, the hydraulic 
properties of the cylinder and servovalve are represented. Many 
system components are represented by nonlinear models; exam- 
ples include the signal quantization of digital to analog convert- 
ers (DACs), pressure versus flow characteristics of valves, and 
potential saturation of the controller command signals. The entire 
system model is shown in Figure 2, and will be referred to in the 
discussion of the modeling and control system design. 

STRUCTURE 

Figure 3 shows the structural model of the pitching mechanism. 
The model was simplified by assuming the cross-member 
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spring constants were large enough to be rigid and, therefore, 
neglected. The structural model can then be simplified to a single 
energy transmission path. Reflecting the parameters to the prime 
mover (cylinder), to eliminate coupling of the individual struc- 
tures, will allow us to solve for the key transfer functions 
defining the equations of motion: 

Cylinder Velocity 
Force Input 

= *± 

F 


FQb) _ A 

P (lb/in. 2 ) 


where A = Area of cvlinder (in. 2 ) 

There is leakage in any hydraulic system and in this model it was 
estimated at 3% of the fluid flow. The leakage will alter the 
damping characteristics for the system. 


leakage 


(3%)(Q in. 3 /sec) 
Available pressure 


' [S 2 + 2(0.05)(1 5.5* lO^S + (15.5M0 3 ) 2 ] 1 
_ 42.8 \ x [S 2 + 2(0.05)(4140)S + (4140) 2 ] J 

S / [S 2 + 2(0.03)(1 5.7* lO^S + (15.7M0 3 ) 2 ] \ 

\ x [S 2 + 2(0.03)(5680)S + (5680) 2 ] j 

Window Position 
Cylinder Position 

-9k 

"Xi 

_ 1.45*10 15 1 

s [ [S 2 + 2(0.05)(15.5* 10 3 )S + (15.5* 10 3 ) 2 ] 

L x [S 2 + 2(0.05)(4140)S + (4140) 2 ] 


where Q = Average fluid rate = Xi*A 

X\ = Cylinder Velocity (188 in./sec) 

A = Piston Area (2.36 in. 2 ) 

Additionally, there is back-velocity fluid flow which is the fluid 
expended in cylinder movement. The back-velocity loop sub- 
tracts from the available fluid flow, much like the back-EMF 
loop in a motor limits the current flowing through it. The back- 
velocity fluid flow is proportional to die speed of the cylinder. 
To drive the cylinder with suitable dynamic response, a high- 
performance servovalve must control the flow of fluid through 
the cylinder. The servovalve controller, internal servovalve 
feedback loops, and mechanical constraints of the servovalve 
must be modeled. The servovalve is modeled with a small- 
signal, second-order transfer function: 

OKs)--, ^"> 2 

S 2 + 2(z)(W„)S + (Wn) 2 


Window Position 
Force Input 

_Ol 

F 

_ 61.9*10 15 1 

S 2 [ [S 2 + 2(0.03)(15.7*10 3 )S + (15.7M0 3 ) 2 ] \ 
\ x [S 2 + 2(0.03)(5680)S + (5680) 2 ] J 

Notice the poles in Ol/Xi cancel with the zeros in Xi/F. This 
property is important in the design of the control system, since it 
enables closing an inner loop (velocity) with zeros in the struc- 
ture thereby aiding in stability. The total combined structure 
transfer function is defined as Ol/F. 


where Kg = Power Spool Flow Gain ((in. 3 /sec)/in.) 
W n = Frequency of servovalve (534 rad/sec) 
z = Servovalve damping coefficient (0.70) 


[S 2 + 2(0.7)(534)S + (534) 2 ] 

Servovalve damping (z) was estimated. To assure 0.70 damping 
is achievable, a bypass valve has been put in parallel with the 
cylinder to add leakage, as necessary, to increase servovalve 
damping. The servovalve has an internal position feedback loop 
closed with a linear velocity differential transformer (LVDT) on 
the first-stage valve spool. This is modeled as an equivalent 
gain, Klvdt. while the servovalve amplifier is modeled as a 
gain, K a (3). 


HYDRAULIC SYSTEM 

The cylinder model and fluid compressibility can now be incor- 
porated into the system model. In the hydraulic system, all fluid 
lines must be rigid and line lengths must be minimized to reduce 
compressibility and system compliance, as well as to account for 
entrained air in the hydraulic fluid. 

P (lb/in. 2 ) _ 4(B) 1 
Q(in. 3 /sec) V T S 

where B = Bulk Fluid Modulus (250,000 psi) 

Vt = Total fluid compressed (27.8 in. 3 ) 

S = LaPlace operator 

The cylinder converts hydraulic fluid flow into pounds of force: 


Klvdt = keW 

where Kq ■ LVDT Demodulator gain (Vdc/Vms) 

Kx = LVDT Spool gain (Vrms/in.) 

Finally, the first-stage valve spool has limited throw which will 
result in system saturation if stroke limits are reached. An 
important nonlinear effect of the servovalve/cylinder system, 
regeneration, must be modeled in a high-response system. The 
cylinder is an energy storage device, much like an electrical 
inductor is a current storage device. The energy of the cylinder is 
stored as momentum in both the piston and the hydraulic fluid. 
The flow through the cylinder cannot be changed instanta- 
neously, and if an attempt is made, a pressure surge will ensue 
much like the voltage surge generated when open-circuiting an 
active inductor. A plot of the flow characteristics of a 
servovalve/cylinder system is shown in Figure 4, where 
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Pl = load (cylinder) pressure, P s = supply pressure (3000 psi), 
Q v = servovalve output flow, and x = first-stage valve 
spool position (3). 


junction in the analog servovalve controller. A gain associated 
with the analog controller is shown after the power amplifier 
saturation block. 


Four quadrants represent both positive and negative pressure, 
Pl, and positive and negative flow, Qy. Quadrants 1 and 3 rep- 
resent the “normal operation” regions where positive pressure 
creates positive flow, and negative pressure creates negative 
flow. Quadrants 2 and 4 represent the regeneration areas. In 
these quadrants the pressure may have reversed polarity, but the 
flow has not yet reversed. This situation occurs when the servo- 
valve receives a command to rapidly change the flow direction 
through the cylinder. As the spool moves and rapidly halts flow 
in one direction, the pressure potential builds up until the cylin- 
der fluid flow changes directions and enters back into quad- 
rants 1 and 3. 

At any flow rate there can be two different pressure values. The 
two equations shown below are used to approximate the previ- 
ously discussed pressure/flow characteristics 

Q>0; Qv = K 

QSO; Qv = K^Ip 

We, therefore, must have a switch in the model governed by the 
direction of fluid flow through the cylinder. This switch will 
choose between the two equations defining the regeneration 
characteristic. Note that the value received from the equations is 
normalized about full flow. For example, if the Pl versus Q v 
calculated flow is 75%, the command flow in the model will be 
multiplied (attenuated) by 0.75 to reflect the pressure drop 
across the servovalve. 


The final system models are the feedback elements consisting of 
a velocity transducer and a position transducer. The velocity 
transducer is a linear velocity transducer (LVT). The motion of a 
permanent magnet induces a voltage proportional to speed in the 
windings surrounding the magnet. The mathematical model is 
simply a gain. The LVT is mounted on the cylinder piston to 
take advantage of the structure zeros enclosed in the velocity 
feedback loop. The position transducer is an incremental 
encoder. The encoder converts angular position of the glass 
disks into digital pulses which are converted into counts in the 
digital position controller. Modeling of the encoder is similar to 
the DA C. A gain is associated with the encoder as angular posi- 
tion in radians is converted into counts in the digital domain. 
Because the analog position is being converted into the digital 
domain, a quantization block is also needed. 

CONTROL ANALYSIS AND DESIGN 

To begin the theoretical control system analysis and design, we 
must tune the innermost loop with all other loops opened. We 
can then close and tune each loop concentrically, until eventually 
all of the loops have been tuned. It is important to note that the 
tuning and analysis will be done with only the linear model 
blocks in place. By eliminating the nonlinear blocks, familiar, 
classical control techniques can be employed. After the optimal 
linear system has been achieved, the nonlinear blocks will be 
added in. These nonlinearities may cause problems with the 
“optimally tuned” linear system and further adjustments would 
then have to be performed. 

Open-Loop Circuit (Xj/Q): 


CONTROLLERS 

A digital controller is used to control the pitch mechanism posi- 
tion loop, while an analog servovalve controller is used for the 
velocity loop. Because the desired motion of the airfoil is a 
ramp, the digital controller must have a ramp generator. The 
controller must also be able to calculate the difference between 
the commanded position and the actual position, which requires 
an equivalent summing junction in software. The controller has 
the capability of compensation in the form of a lead and lag fil- 
ter. This capability enables placing a zero and/or pole into the 
position loop, as well as an integrator to help remove steady- 
state errors. 

Because of the nature of digital control and the sampling of the 
analog system, there is an inherent time delay associated with the 
controller’s operation. If the sampling rate is not suitably high, 
the time delay could be detrimental to the stability of the system 
and so it is modeled. The system is analog, so the digital con- 
troller must convert its output to the analog domain. A DAC 
performs the conversion which creates a quantized analog sig- 
nal. The quantization could potentially excite a mode in the sys- 
tem and is important to model. The DAC also has a gain associ- 
ated with it The gain is a coupling ratio of counts in the digital 
domain to voltage in the analog world (6). 

The output voltage from the DAC is fed to an analog voltage 
power amplifier which amplifies the signal sent to the servo 
controller. Because this system requires high dynamic response, 
it is possible that the power amplifier will saturate if the input 
signal exceeds ±10 V which must be included in the model. The 
analog command voltage is then fed into the velocity summing 


This discussion begins with the cylinder, area, structure, and 
leakage feedback blocks of the system model. The theoretical 
poles and zeros of this system are shown below: 


Zeros 

Comolex Roots 

Frequency, 

rad/sec 

Damping, 

_z 

Structure: 

-777 3 15,520 

15,520 

0.05 

Structure: 

-207 34,140 

4,140 

0.05 


Frequency, Damping, 

Eoks Complex Roots -tad/sec z 


Structure: 

-472 315,730 

15,730 

0.03 

Structure: 

-170 35,670 

5,670 

0.03 

Cylinder 

-425 

425 

1.0 

Velocity 

integrator: 

0 

0 

1.0 


The structure poles are lightly damped, although at a high natural 
frequency, while the cylinder has a low natural frequency. To 
get the desired performance from the system, the cylinder fre- 
quency will have to be improved. The open-loop Bode plot of 
this system is shown in Figure 5. 


Back- Velocity Loop (Xi/Q): 

The next loop to close is the back-velocity loop which feeds 
back the cylinder speed. The Bode plot is shown in Figure 6. 
The new system roots are shown below: 
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Zeros 

Complex Roots 

Frequency, 

rad/sec 

Damping, 

z 

Structure: 

-777 ±) 15,520 

15,520 

0.05 

Structure: 

-207 ±j4,140 

4,140 

0.05 

Poles 

Structure: 

Complex Roots 

Frequency, 

rad/sec 

Damping, 
. z 

-460 ij!5,740 

15,750 

0.03 

Structure: 

-174 ±j6,050 

6,050 

0.03 

Cylinder/ 

integrator 

-220 ±j 1,960 

1,970 

0.11 


The velocity integrator, which came from the Xi/F transfer 
function, combined with the cylinder pole and became a complex 
root pair at a higher natural frequency, which is achieved at the 
expense of cylinder damping. The lower frequency structure 
pole moved out to a slightly higher frequency, and the high- 
frequency pole did not move. 

The next step is to add the servovalve with its associated LVDT 
loop and amplifier to the previous system. The open-loop Bode 
plot is shown in Figure 7. The resultant system frequencies are 
shown below: 


Zeros 

Complex Roots 

Structure: 

-777 $ 15,520 

Structure: 

-207 $4,140 

Poles 

ComplexRoots 

Structure: 

-460 $ 15,740 

Structure: 

-174 ±j6,050 

Cylinder/ 

integrator 

-220 $ 1,960 

Servovalve: 

-377 $378 


Frequency, 

rad/sec 

Damping, 

z 

15,520 

0.05 

4,140 

0.05 

Frequency, 

Damping, 

rad/sec 

z 

15,750 

0.03 

6,050 

0.03 

1,970 

0.11 

534 

0.70 


Velocity-Feedback Loop (SCi/V): 

The velocity loop will force the cylinder to accurately track the 
velocity command. To appropriately choose the velocity loop 
gain, a root locus plot is created to track the poles’ locations as a 
function of velocity loop gain. The loop gain must be increased 
in the velocity loop to increase the servovalve pole frequency 
and get the necessary performance from the servovalve. Notice, 
however, in the previous Bode plot, that the phase and gain 
margins are negative, indicating an unstable system with unity 
feedback in the velocity loop. The velocity-feedback transducer 
gain must, therefore, be less than unity. The root locus in Fig- 
ure 8 confirms this assertion, showing the cylinder poles head- 
ing toward the right-hand plane as velocity-feedback gain 
approaches unity. 

A feedback gain for the velocity transducer was chosen and the 
resultant roots arc shown below: 

Frequency, Damping, 
Zeros Complex Roots rad/scc 2 


Poles 

Complex Roots 

Frequency, 

rad/sec_ 

Damping, 

z 

Structure: 

-460 $ 1 5,740 

15,750 

0.03 

Structure: 

-174 $6,050 

6,050 

0.03 

Cylinder/ 

integrator 

-212 $1,920 

1,930 

0.11 

Servovalve: 

-386 $587 

703 

0.55 


Notice that the servovalve frequency was improved at the 
expense of damping. The Bode plot is shown in Figure 9. As 
stated previously, the velocity and position transducers are 
located in two different places. The difference in location of the 
transducers can be explained by looking at the transfer function, 
Xi/F. This transfer function has two zero pairs in the numerator 
of the transfer function. These zeros, enclosed in the velocity 
loop, tend to attract system poles, keeping the poles from enter- 
ing the right-hand plane, and thereby becoming unstable. 
Because these poles are restrained, the velocity loop gain can be 
increased to get the velocity of the device to more accurately 
track the velocity command. If the velocity transducer had been 
put at the airfoil instead of the cylinder, the zeros would have 
been cancelled by the poles in the 0\JX\ transfer function. The 
system poles would not then have been restrained to the left- 
hand plane. The velocity transducer is, therefore, located on the 
cylinder in the velocity-feedback loop, while the position loop 
transducer is located at the actual load, to assure that the airfoil 
arrives at its final position as desired. 

Position Feedback Loop (Ol/Q: 

The lead and lag filters in the position loop of the digital con- 
troller are initially set to zero to best tune this system prior to 
compensation. The position loop gain was chosen to achieve the 
best dynamic response. In this case the servovalve poles are 
heading toward the right-hand plane, limiting the amount of 
position loop gain. Resultant roots are shown below: 


Zeros 

Comolex Roots 

Frequency, 

rad/sec 

Damping, 

z 

Structure: 

-777 $15,520 

15,520 

0.05 

Structure: 

-207 $4,140 

4,140 

0.05 

Poles 

Complex Roots 

Frequency, 

rad/sec 

Damping, 

z 

Structure: 

-460 $15,740 

15,750 

0.03 

Structure: 

-174 $6,050 

6,050 

0.03 

Cylinder/ 

integrator 

-247 $1,920 

1,930 

0.13 

Servovalve: 

-80.2 $608 

613 

0.13 

Position 

integrator 

-544 

544 

1.0 


It can be seen that the servovalve poles migrated toward the 
right-hand plane of the root locus, thereby decreasing the fre- 
quency and damping of that pole pair. The real pole is also at a 
low frequency, resulting in decreased system bandwidth. The 
Bode plot for the uncompensated position loop is shown in 
Figure 10. 

Compensation (Ol/Q: 


Structure: -777 ±j 15,520 15,520 

Structure: -207 ±j4,140 4,140 


0.05 The position controller software lead compensation can now be 

0.05 used to improve the damping of the servovalve. The zero is 

placed at 900 rad/sec, while the accompanying compensation 
pole is placed at 2000 rad/sec. The compensation integrator in 
the software will not be used. Although the integrator helps to 
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reduce steady-state error in the system, it also introduces a 90° 
phase lag which decreases phase margin and, therefore, system 
damping. The resultant root locus plot is shown in Figure 11. 
The system roots after compensation are shown below: 


Zeros 

Complex Roots 

Frequency, 

rad/sec 

Damping, 

z 

Structure: 

-776 ±j 15, 520 

15,520 

0.05 

Structure: 

-207 ±j4,l40 

4,140 

0.05 

Lead 

compensator: 

-900 

900 

1.0 

Poles 

Complex Roots 

Frequency, 

rad/sec 

Damping, 

z 

Structure: 

-460 ±j 15,740 

15,750 

0.03 

Structure: 

-173 ±j6,050 

6,050 

0.03 

Cylinder/ 

integrator 

-248 ±j 1,900 

1,910 

0.13 

Servovalve: 

-189±j639 

666 

0.28 

Position 

integrator: 

-297 

297 

1.0 

Lag 

compensator: 

-2.030 

2,030 

1.0 


Improvement in the servovalve frequency and damping is 
noticeable. The servovalve natural frequency has been improved 
from 613 rad/sec to 666 rad/sec (10%), and the damping 
improved from 0.13 to 0.28 (115%). The Bode plot for the 
compensated position loop is shown in Figure 12. 

Nonlinear Elements: 

Now that the linear control system model has been tuned, the 
nonlinear model blocks can be inserted and the actual predicted 
results viewed. Various system responses are displayed in Fig- 
ure 13. An important plot to note is the position controller 
amplifier where the output goes into saturation. When the posi- 
tion loop is saturated, the effect is to open the position loop 
feedback, giving the velocity-feedback loop control of the sys- 
tem. The response of the system is therefore determined by the 
poles in the velocity loop while the position loop is saturated. 

The effect of this amplifier saturation can be seen in the plot of 
flow response. Upon close inspection, it can be seen that the 
frequency of the ringing at high-flow command is different from 
the ringing frequency at zero flow command. The frequency of 
oscillations during high flow, while the position controller is 
saturated, is 105 Hz; while at zero flow command the oscillation 
frequency is only 95 Hz. These ringing frequencies are deter- 
mined by the natural frequency of the dominant pole(s) in the 
velocity and position loops, respectively. The servovalve is the 
dominant pole in both the velocity loop and the position loop. 
The output ringing frequency found during position controller 
saturation corresponds to the servovalve frequency in the veloc- 
ity loop. The output ringing frequency found while the position 
controller was not saturated corresponds to the servovalve fre- 
quency in the position loop. 


Note that in the pressure response plot, the pressure during 
acceleration reaches 900 psi, which is less than the deceleration 
pressure of 1300 psi. This difference demonstrates the back- 
velocity loop’s ability to aid in stopping the cylinder motion. 
During the acceleration, the back-velocity loop is fighting the 
motion all the way, impeding the cylinder’s ability to accelerate. 
During deceleration, the back-velocity loop is still trying to 
impede the motion of the cylinder, aiding in slowing the cylinder 
to a stop. 

The high-frequency shaking occurring in the cylinder accelera- 
tion plot (X) appears to be one of the structure poles. The 
shaking frequency is approximately 1000 Hz which matches the 
lower frequency structural pole of 962 Hz. This plot indicates 
the potential for exciting some of the structure poles. However, 
in the position-response plot the shaking is integrated twice and 
dramatically reduced. 

CONCLUSION 

A system model encompassing a wide range of real-world 
effects showed the ultimate system design parameters of 60° 
pitch angle at 3600°/sec pitch rate can be theoretically achieved. 
The acceleration of the airfoil was within 0.006 sec and 10° of 
pitch angle. Although the required design performance was 
demonstrated, there is little room for improvement in the system 
response. All of the system variables were maximized for 
linearity of pitch rate, minimum overshoot, and appropriate 
acceleration rates. 
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Figure 1 - Constant pitch rate apparatus. 
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Figure 2 - Control system model. 
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Figure 3 - Mechanical structure model. 


Figure 4 - Servovalve pressure-flow characteristics. 
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Figure 5 - Bode plot of open loop structure (Xi/Q). 
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Figure 6 - Bode plot of back-velocity loop (Xj/Q). 
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Figure 9 - Bode plot of velocity-feedback loop (Xj /V). 
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Figure 7 - Bode plot of structure and servovalve. 
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Figure 10 - Bode plot of uncompensated position loop. 



Figure 8 - Locus of roots for velocity feedback. 
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